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^ ', 1. Introduction 

ft . 

The electron temperature gradient (ETG) driven mode [U [2j [3] is often considered as 
a possible candidate for the explanation of electron thermal transport through internal 
transport barriers when ion temperature gradient (ITG) turbulent fluctuations are 
suppressed by E x B shear flow. At the same time, small scale ETG fluctuations have 
the typical spatial scale of the order of electron gyroradius p e , and are less susceptible 
to quenching by shearing [U G3 EJ- It is widely thought that drift-wave-type turbulence 
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can excite zonal flows which are associated with azimuthally symmetric band-like shear 
flows that depend only on the radial coordinate IB] • Zonal flows play a crucial role in 
regulating drift-wave turbulence and transport in tokamaks. It is now quite clear that 
zonal flows are generated by modulational instability of drift waves P, QUI EEH 021 H3] . 
For ETG driven turbulence, the excitation of zonal flow was considered in [T3] (the case 
of broad turbulent spectrum) and [15] (the four- wave coupling scheme). 

It is well known that the presence of shear flow give rise not only to instability 
of the sheared layer (Kelvin-Helmholtz instability), but also to stabilization of other 
instabilities (ITG driven modes, resistive interchange modes etc.) [3 HI E3, EHl EH] . Up 
to now, the stabilizing effect on the linear instabilities has been considered only for the 
case of mean smooth flows. Note, that zonal flow shearing differs from that for mean 
flow shearing on account of the complexity of the flow pattern. In contrast to smooth, 
static mean flows, the zonal flow patterns can be expected to have finite correlation time 
and complex, possibly random, spatial structure [71 120]. 
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In the present work we consider the influence of zonal flow on the linear stability 
of ETG drift modes. We show that the presence of random shearing by zonal flows 
strongly affects the linear stability of ETG modes and has a stabilizing effect. If the 
mean square amplitude of zonal flow exceeds some critical value, the linear instability 
of ETG modes is suppressed for all poloidal wave numbers k y . 



2. Basic equations 

Assuming a slab two-dimensional geometry charge quasineutrality and the adiabatic ion 
responce, we consider the following simplified model describing curvature driven ETG 
turbulence and including viscosity and thermal diffusivity [5] 

^ + v E - VJ (<p-A ± <p) + —{<p + p) + uA 2 <p = 0, (1) 
d \ d 

^ + v E - V )p-r— y;-xAp = 0, (2) 

where vj = [z x V<£>] is E x B drift velocity, ip and P are the normalized electrostatic 
potential and plasma pressure respectively, and {A, B} = d x Ad y B — d y Ad x B is the 
Jakobian. The variables in equations (CQ) and ([2j) have been rescaled as follows 
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t — e*iUJBit', 

x', y' and t' being the original physical coordinates (with x' the poloidal and y' the 
radial coordinate), 

_ p s ^/t _ p s y/r _ p s ^/t 

J-'n ^B L'p 

where L n , Lb and L p are the background gradient scales for the density, magnetic field 
and pressure respectively, p s is the ion gyroradius calculated at the electron temperature 
T e , and r = Tj/T e . The effect of zonal flow can be included by assuming the E x B drift 
velocity of 

v E = v(x)y + [z x Vtp], (3) 

where the first term accounts for zonal flow. Neglecting nonlinear terms in equations 
([I]) and (j2]), taking into account equation (j3J) and representing 

(p(x, y, t) = exp(ik y y — hot), 



p(x, y, t) = P(x) exp(ik y y — hot), 
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one can obtain 

(-iu + ik y v)(l + k 2 y - + iky(* + P) + viky + ^ - 2*^)* = 0, (4) 

(-iw + - i V $ + x(-k 2 y + ~[^) p = °- (5) 

In the inviscid limit and absence of zonal flow, after taking <&(x), P(x) ~ exp(ik x x), 
equations (jlj) and ([SD give the dispersion relation for ETG modes 

k,i 



2(F + 1)l l±Vl-4r(^ + l)j, (6) 

where k 2 = k 2 + k 2 , and the plus sign describes the drift waves dispersion, while the 
minus sign corresponds to the dispersion of the convective cells. Equation ([6]) predicts 
instability with the growth rate 

kn 
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2(k 2 + l 
if 4r(k 2 + 1) > 1. 



V4r(P + 1) - 1, (7) 



3. ETG drift modes in the presence of zonal flow 

We assume that the profile of zonal flow v (x) is the random function and v(x) is a zero 
mean, (v(x)) = 0, real homogeneous Gaussian process with correlation function 

(v(x)v(x')) = D(x-x'), (8) 

where (...) means statistical averaging. It is assumed, that the intensity of the noise is 
small, D(x) <C 1. When neglecting v 2 term, and in the inviscid limit, equations (jl]) and 
([5]) reduce to the equation 

w2 (^-^ 1 ) $ + W f-^(^"^ 1 )} $ "* + r)$ = ' (9) 
which after the Fourier transforming (d/dx —>■ —iq) can be rewritten as 

G \p)$(p) + J k y [k y - 2u{q 2 2 + k 2 y + l))v{ qi )<f>{u, q 2 , k y )S(q -q x - q 2 )d qi dq 2 = 0, (10) 

where p = (u,q,k y ) and 

G \p) = u 2 {q 2 + k 2 y + 1) - u;k y + k 2 y r. (11) 
To avoid the cumbersome expressions, we write equation fflOl) in the symbolic form 

G^ipMp) + l{p,Pi,P2)v{pi)$(p 2 ) = 0, (12) 

where 

1f(p,Pi,P2) = x{vi,V2)5{p - P\ - Pi), (13) 

and integration over repeated indices is assumed, as usual. Introducing a source rj(p) in 
the right hand side of equation (1121) and taking the functional derivative 5/5rj(p'), we 
have 

G 1 (p)G(p,p') +i(p,Pi,p 2 )v( Pl )G(p 2 ,p') = S(p-p'), (14) 
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where G(p,p') = 5<&(p) / 5rj(p') is the Green function, and we have taken into account 
that Srj(p)/6r](p') = 5(p — p'). The function G (p,p') is a free (i. e. in the absence of 
the noise v) Green function. Representing the Green function as a sum of the average 
and fluctuating parts 

G(p,p') = (G(p,p , ))+G(p,p'), (15) 
substituting it into equation ( |T4l) and averaging, one can obtain 

G 1 (p)(G(p,p'))+ 1 (p,p 1 ,p 2 )(v( Pl )G(p 2 ,p')) = S(p-p'), (16) 
Substracting equation ffTB"]) from equation ffT4"l) . we get 

G 1 (p)G(p 1 p') + 7(P,PhP2)v(pi)(G(p 2 ,p')) 

+l(p,PuP2)[v(pi)G(p 2 ,p') - {v{px)G{p2iP'))\ = 0. (17) 

In the Bourrett approximation |21j . which is justified when the intensity of the noise is 
small enough, we can neglect the term v G — (vG) and get for the fluctuating part of the 
Green function 

G(p,p') = -G {p) 1 (p,p 1 ,p 2 )v( Pl )(G(p 2 ,p')}. (18) 
Inserting this expression into equation ( jl~6l) . we obtain 

G 1 (p)(G(p,p')} -i(p,PuP2h(p' 1 ,p' 2 ,P2)Go(p2)(v(pi)v(p , 1 )}(G(p' 2 ,p')} = 5{p-p'). (19) 

As follows from equation ([8]), in the wave number domain the correlator has the form 
(v{p)v{p r )) = I{p)8{p + p'). Due to homogeneity of the random process, the Green 
function in equation (|T9j) has the structure (G(p,p')) = G{p)5{p — p'). Making use of 
this and equation (JT3l) . one can perform some integrations in equation (fl~9l) and finally 

we get 

GP Gq 1 (p) - x(p-px,p 1 )x(p 1 -p,p)G (pi)I(p-pi)' ^ 
The poles of the Green function G(p) determine the spectrum of elementary excitations 
and the corresponding dispersion relation is 

G 1 (p) - x(p- Pi,Pi)x{Pi -P,p)G (px)I(p-p 1 ) = 0. (21) 

In the absence of zonal flow we have Gq 1 (p) = and recover the previous result (J6]) 
with k x = q. In what follows we consider the case when the profile of the zonal flow is 
described by the random function v{x) which has the form 

v (x) = Vq cos(g £ + (22) 

where the random amplitude v o is a zero mean, normally distributed value with variance 
a 2 , and the random phase $ is uniformly distributed between and 2n. The correlation 
function (JHJ) of such a process is D(x) = (c 2 /2) cos(go^) or, in the wave number domain 

2 

I{q) = ^[S(q-qo) + S(q + q )]. (23) 

In this case the noise has an infinite correlation length and is concentrated at the wave 
number q so that the characteristic scale length of the zonal flow is of the order ~ q^ 1 . 
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Figure 1. (A) The instability growth rate as a function of the root mean square zonal 
flow amplitude a (for different values of the poloidal wave number k y ). (B) The growth 
rate as a function of the poloidal wave number k y (for different values of a). 



As follows from eqution ([7]), in the absence of zonal flow the maximum growth rate for 
the fixed k y is achieved at k x = 0. Thus, to simplify calculations, we will consider the 
influence of zonal flow on the most unstable modes and put q = 0. Then, substituting 
equation (1231) into equation (1211) one can obtain the dispersion relation 

cito 4 — C2UJ 3 + C3U 2 — C40J + C5 = 0, (24) 



where 



Cl = (^ + l)(g 2 + ^ + l), (25) 

c 2 = k y (q 2 + 2k 2 y + 2), (26) 

C3 = k 2 y (q 2 + k 2 y + l)[r- 2a\k 2 y + 1)] + k 2 y [l + r(k 2 y + 1)], (27) 

Ci = 2k 3 y r-ky(q 2 + 2k 2 y + 2), (28) 



a 



Cb = A{ ( t ~ 2 - - ) • (29) 



Figure Ufa) shows the growth rate 7 of the ETG driven mode as a function of the root 
mean square zonal flow amplitude a for different values of the poloidal wave number 
k y . The characteristic wave number of zonal flow used in the calculations is c/o = 0.1. 
The parameter r has been fixed at r = 0.25 so that, as follows from equation ([7j), in the 
absence of random shearing there is a linear instability for all k y 7^ 0. It is seen that 
the growth rate decreases with increasing a and vanishes above some value of a which 
depends on k y . Thus, zonal flow shearing stabilizes the instability of ETG modes. In 
figure Mjo) we plot the dependence of the growth rate on the poloidal wave number k y 
for different values of the root mean square zonal flow amplitude a. In the absence of 
zonal flow (er = 0) the growth rate 7 increases with increasing k y and saturates at the 
level ■jmax = y/r- The presence of shearing (er 7^ 0) changes the situation drastically. 
The reduced growth rate initially increases as a function of k y from some k Vy i and then 
decreases, becoming zero at some k y ^- For not too large er, the linear instability is 
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Figure 2. Maximum growth rate versus k y for different values of q in equation (|30|) 
and (A) zonal flow amplitude vq — 0.5; (B) zonal flow amplitude vq = 1. 



restricted to the region k y> i < k y < k yj2 - The lower k y> i and upper k y ^ boundary 
values increase and decrease respectively as the root mean square zonal flow amplitude 
a increases and above some critical value a cr ~ 0.36 (for r = 0.25) the instability of 
ETG modes is suppressed for all poloidal wave numbers k y . An estimate for the critical 
value a cr can be obtained from equation (12^1) . The growth rate remains zero as k y — > oo 
(the most dangerous case) if a = a cr . From equation (12^1) one can see that u scales as 
k^ 1 as k y — ^ oo. Then, we can get the estimate a cr = y/2r. This theoretical prediction 
for the dependence of the critical value of root mean square zonal flow amplitude on r 
is in very good agreement with numerical results. 

Next, we consider the case when the zonal flow profile is deterministic and has the 
form v(x) = focos(gx). In addition, we include the effects of viscosity and thermal 
diffusivity. Then, equations Q and (jSJ) can be rewritten as an eigenvalue problem 





where 



k y A- x B 



d 2 v - 

A = l + k 2 y -— , B = 1 + v cos(qx)A-i— D, 

ijjjb y 



(30) 



d 2 d A x d 2 

w + ^ ' c = Vo cos ( qx) ~ lxky + l Td^- 



"y 

Employing a finite differencing approximation, we numerically solved the eigenvalue 
problem (|30|) . The dissipative coefficients have been fixed at v = x = 0.01. The 
instability growth rate is plotted in figure [2] for two values of the zonal flow amplitude 
Vq and different values of q (for r = 0.5). It is seen that the presence of zonal flow 
reduces the growth rate though the stabilizing effect manifests itself not so sharply as 
in the case of random shearing. 
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4. Conclusion 



In conclusion, we have investigated the influence of zonal flows on the linear instability of 
ETG driven modes. Random and deterministic cos - like profiles of the zonal flow have 
been considered. For the random profile of zonal flow, we have obtained in the Bourrett 
approximation the dispersion relation for ETG modes in the presence of shearing. We 
have shown that the presence of random shearing caused by zonal flow has a strong 
stabilizing effect on the ETG driven mode destabilized by the temperature and pressure 
gradients. If the mean square amplitude of zonal flow exceeds some critical value, the 
linear instability of ETG modes is suppressed for all poloidal wave numbers k y . 
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